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Abstract 

We propose a first-principle method for evaluations of the time-dependent electron distribution 
function of excited electrons in the conduction band of semiconductors. The method takes into ac- 
count the excitations of electrons by external source and the relaxation to the bottom of conduction 
band via electron-phonon coupling. The methods permits calculations of the non-equilibrium elec- 
tron distribution function, the quasi-stationary distribution function with steady-in-time source 
of light, the time of setting of the quasi-stationary distribution and the time of energy loss via 
relaxation to the bottom of conduction band. The actual calculations have been performed for 
titanium dioxide in the anatase structure and zinc oxide in the wurtzite structure. We find that 
the quasi-stationary electron distribution function for ZnO is a fermi-like curve that rises linearly 
with increasing excitation energy whereas the analogous curve for anatase consists of a main peak 
and a shoulder. The calculations demonstrate that the relaxation of excited electrons and the set- 
ting of the quasi-stationary distribution occur within the time no more than 500 fsec for ZnO and 
100 fsec for anatase. We also discuss the applicability of the effective phonon model with energy- 
independent electron-phonon transition probability. We find that the model only reproduces the 
trends in changing of the characteristic times whereas the precision of such calculations is not high. 
The rate of energy transfer to phonons at the quasi-stationary electron distribution also have been 
evaluated and the effect of this transfer on the photocatalyses has been discussed. We found that 
for ZnO this rate is about 5 times less than in anatase. 



I. INTRODUCTION 

The zinc oxide in wurtzite structure and the titanium dioxide in anatase structure are the 
semiconductors of great interest from many viewpoints. Both demonstrate photocatalytic 
activity in the UV region of the sunhght spectra and are widely used as basic compounds 
in attempts to create the photocatalysts active in the visual sunlight that can be used for 
cleaning the environment from organic pollutants and pathogenic bacteria^"^. They are 
also perspective materials for applications in solar energetics and random access memory 
devices^"-. The photocatalytic activity, magnetic properties and charge transport in these 
oxides to a great extent depend on the dynamics of the excited states, therefore in many 
works the characteristics of fast electron and hole dynamics have been studied. In the 
worka^i"— the relaxation rates of excited electrons in ZnO via electron-phonon coupling have 
been estimated whereas in the papers^'— the characteristics of relaxation rates for anatase 
have been evaluated. 

The comprehension of the results obtained in the experimental works requires an applica- 
tion of some theoretical models that should incorporate both the conditions of experiments 
and the characteristics of electron dynamics determined by the properties of the solids un- 
der study. Such models have been successfully developed and widely applied for the excited 
electron relaxation in metals^^"— . In the case of semiconductors this aim is, however, far 
from being achieved, partly because only recently the rigorous methods for evaluations of 
the characteristics of electron dynamics have been developed. The main phenomenon re- 
sponsible for the relaxation of the low-energy excited electrons in semi-conductors (whose 
energy with respect to the bottom of the conduction band is no more than the band gap) is 
the electron-phonon coupling. A number of first-principle method for the calculations of the 
parameters of this coupling have been developed that permit to perform subsequent calcu- 
lations for the temporal characteristics of electron relaxation. Possibly, the most promising 
of them is provided by the density-functional perturbation theory^"^ implemented in the 
pseudo-potential Quantum Espresso (QE) computer code. Recently this approach has been 
applied to the evaluations of the characteristics of electron relaxation in ZnO^^, anatase and 
rutile^^, germaniunt^i^ . The calculations have demonstrated a reasonable correspondence to 
experimental data and helped to understand some fine relations between the band structure 
and electron dynamics in the compounds. 



The shortage of the evaluations of the Refs.— 12^ was that they had been performed for 
a case of a single excited electron in the conduction band that falls down to the bottom of 
the band loosing the energy via the electron-phonon coupling. More rigorous estimations 
of the electron dynamics should include also the emergence of new electrons excited by the 
external light source and also the effect of filling the energy level under consideration with 
the electrons falling from the higher levels. The method of such evaluations is proposed 
in the current paper. Although the concrete results are obtained for ZnO and Ti02, the 
proposed method does not have features preventing its application to different kinds of semi- 
conductors. The method is based on the fully first-principle calculations of the electron- 
phonon coupling constants. So it does not have simplifying assumptions on the energy 
dependence of the rate of electron-phonon relaxation. It is worth noting that the neglect of 
the energy dependence of the electron-phonon interaction is a conventional approximation, 
as it is e.g. with the Frohlich's hamiltonian^^"— . So, together with the characteristics 
of electron relaxation in ZnO and Ti02, we discuss also the applicability of the energy- 
independent electron-phonon coupling and find that in the current case it is of very limited 
accuracy. 

II. METHOD OF CALCULATIONS 

A. Quasi-stationary electron distribution 

We start with the following Boltzmann equation for the time-dependent distribution 
function /(t, e) of excited electrons (hereafter EDF) at the energy e above the bottom of the 
conduction band of a semiconductor (hereafter - excess energy) 

df(t e) Z"'^™ 

•'\' ' = (1 - nit, e)) / nit, e + nuj)Fie + hu, u)duj- 
ot Jo 

nit, e) r™ F(e, u)duil - n{t, e - huj)) + ^^'"1*^^'"^ (1) 

JQ at 

In this equation nit,e) = fit,e)/Ni^e), where A^(e) is the density of electronic states, is the 
population of a single band state at the e level. First integral describes, in a momentum- 
averaged manner, the process of coming of electrons to the level e from all the electronic 
states higher by the energy hu, the process accompanied by the emission of phonons whose 
maximum frequency is Um- F{e+huj, u), the spectral function of electron-phonon interaction. 



is the probability for electrons to fall from all the band states with energy e + hu to the 
states at the energy e emitting phonons with energy hu satisfying the conditions of energy 
and momentum conservation. Similarly, the second integral describes the process of electron 
transitions from the e level to the levels lower by the phonon energy hu. We discuss only 
the case of a small intensity of irradiation, so we assume that the n{t, e) values are small 
and omit the factors 1 — n{t,e). The spectral function F{e,u), satisfying to energy and 
momentum conservation, has the form 

F(e, w) = ^ ^ 6{e - enk)P{nk, n'k' , o-q)5k-k'±q5(w - ujaq)6{e - hu - e„/k') (2) 

qkk' nn'a 

where P(nk, n'k.', aq) is the probability of a single electron transition between the electronic 
band states \nk) and \n'k') accompanied by the emission or absorption of the photon with 
frequency Wcrq- The calculations of such spectral function, which for the Fermi level in 
metals is proportional to the well-known Eliashberg function, were extensively discussed in 
literature^i^S. In accordance to the 'golden Fermi rule', this probability, per unit of time, is 
determined as 

P{nk,n'k',crq) = {2n/h)\{nk\AV^^\n'k')\^5{enk - e^'k') (3) 

Here the value AV^g. is the variation of the self-consistent potential in crystal caused by the 
displacement mode of the phonon, and the value (nk|AV^ty|n'k') is the matrix element of the 
electron-phonon interaction. The energy of phonon, small in comparison with the energy of 
electronic states, is neglected here. The ways of calculating this matrix element also have 
been well discussed^^. 

The term <9/j„st(t, e)/9t in Eq. [1] describes the instantaneous electron distribution in the 
conduction band produced by the external source of excited electrons. One can take for this 
term the approximation 

^^^"gf''^ = So{t)Sie) (4) 

where the factor So(t), which in general can be time-dependent, is the concentration of the ex- 
cited electrons determined by the power of the light radiation, and S{e) is the instantaneous 
spectral function ( hereafter lEDF) which describes the probability for an excited electron to 
have the excess energy e. Naturally, lEDF should be normalized to unity: J^ S{uj)du = 1. 
If one linearizes the energy dependence of the F{e + hu, u) with w as a small parameter 
then Eq. [1] is written in the form 



^ = ||iV-.(.)/(M)BWl + ^i=^ (5) 

Here the function 

5(e) = / huF{e,uj)du (6) 

Jo 

is the energy lost by the electrons at the level e via phonon emission. So the value which 
we define as huo = B{t)/N{t) is the averaged energy lost by one electron in the process 
of phonon emission; a first-principle approach to the calculations of this value has been 
demonstrated in Refs.— "^ (where it has been denoted as Ae). 

Eq. [5] describes the time evolution of the f{t,e) EDF when the light source produces 
an instantaneous distribution of excited electrons in the conduction band. We discuss here 
the quasi-stationary case, when the temporal evolution of the external pulse is slow and 
dSo{t)/dt « So{t)B{e). Taking dfit,e)/dt = one has the solution of the Eq. [5] in the 
form 

f{t,e) = Soit)j^J^^"^S{e')de' (7) 

where e^ is the highest excess energy of the excited electrons. 

The lEDF S{uj) can be obtained from the electronic band structure calculations. Namely, 
if the energy of the quantum of optical excitation is Egxc, then for the excess energy e one 
has to sum the probabilities of all direct excitations from the electronic states at the energy 
e — Eexc to the states at the energy e. Hence, the un-normalized lEDF is 

S{e) = "^Sie- enk)T{nk, n'k)6{e - E^xc - e„'k) (8) 

knn' 

where T(nk, n'k) is the probability of the transition between the states |nk) and |n'k). 
(In practical calculations we replace the (5-functions with the normalized to unity gaussians 
whose width at the half-maximum is 0.01 eV.) The transition probability also can be evalu- 
ated basing on the first-order perturbation theory. In order to calculate the matrix elements 
of the Eq. ([8]) we take advantage of the atomic sphere approximation'^^. In this approxima- 
tion the integration over the space of a crystal is replaced with the integration over atomic 
spheres. In every atomic sphere we take for the perturbation the dipole approximation^^. 
With this approximation the hamiltonian of the interaction of an atom with the electric field 
of light Fj{u) has the form 

Hs = eD5E(u;) (9) 



Here eD^ = eJ2j i"f ? where r^ is the radius with respect to the center of the given atom, is 
the operator of the dipole moment of the atom. Hence, the interaction of an electron with 
the field is 

H{iu) = eY,r^eEE{u) (10) 

s 

We consider the case of interaction of light with a polycrystal, so we have for the angle- 
averaged transition rate 



T{nh, n'h) 



(11) 



Here the overlap of atomic spheres is neglected, and the coefficient 1/3 emerges because of 
the averaging over the angle between the vector r'^ and the directional vector of the field eq, 
see details iniS. 

B. Electron-phonon energy loss time and the time of EDF setting 

The value F(e + Tiuj.uj) in Eq. [H according to the definition ([2]), is the probability of 
transitions from all the electronic states at the energy e + hu to all the states at the energy 
e. So we apply to this value the approximation 

F(e + hu, u) = N{e + hio)P{e)N{e)S{io - uo) (12) 

where -P(e) is the momentum-averaged probability of a single transition at the electron 
energy e and phonon energy huQ. The energy dependence of P is often neglected, as it takes 
place in numerous works with Frohlich's electron-phonon interaction, see e.g.—"—. We will 
show in the next section that the energy dependence of -P(e) can be omitted only for the 
aims of interpretations. In fact this dependence is not negligible. Employing approximation 
( IT2II and introducing for the external source definition 

dninst(t,e) _ dfinstJt , e) 
dt ~ dt 

we can rewrite Eq. ([T]) in the form 



-/N{e) (13) 



^^^ = n{t, e + huo)Pie)Nie + huo) - n{t, e)P(e)iV(e - ^c^o)] + ^""^"gf ' '^ (14) 



After linearization near the energy e the equation becomes 

We introduce a new variable 

g(t,e)=iV2(e)n(t,e) (16) 

and get for this value the equation 

^ = huoPie)Nie)^l^ + ^^^#1^ (17) 

dt ^ ^ ^ ^ de dt ^ ^ 

One can show that in the absence of light the formal solution of Eq. (fT7|) at the energy E 
can be written as 

q{t, E) = ${t + [hu,]-^ r de{P{e)N{e))-^ (18) 

JEi 

where $ is an arbitrary function; this function and the energy Ei have to be chosen in order 
to satisfy initial conditions. We assume for the initial conditions that at t = the electron 
is excited by the external source to the level Ei, and afterwards the source is switched off. 
It is easy to check that these conditions are satisfied if $ is a 5-like function, so 

n{t,E) = N-^5{t+{hujQ]-^ j de{P{e)N{e))-^} (19) 

JEi 

The equation 

t + [hujo]-^ r de{P{e)N{t))-^ = (20) 

JEi 

is then that of relaxation of the excited electron. Hence for the rate of the energy relaxation 
we have 

dE/dt = -huoP{E)N{E) (21) 

and the energy-loss time, that is the time necessary for the electron to fall from the level Ei 
to the bottom of the conduction band, is 

TiiE,:) = (huo)-' r\P{e)N{e))-He (22) 

JO 

In real calculations the lower limit of integration has to be replaced with the maximum 
energy of phonons Ttujm since at lower energy the distribution function is determined by 
different mechanisms, mainly by the electron-hole recombination. 
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Eq. (122|) demonstrates the relation between the time of energy loss and the density of 
states. One can reveal also the relation between the energy-loss time in the current paper 
and that in the previous works^i2^. In the cited papers the rate r(e) of the electron-phonon 
relaxation has been defined via 

r(e) = P{e)N{e). (23) 

So the r value is the transition probability from the electronic state at energy e to all 
the states inside the energy interval from e — hUm to e, and the equation for the energy-loss 
time is 

TiiE,) = (htoo)-' rT-\e)de (24) 

Jo 

Since the value Huq = Ae is almost energy-independent^i2^, this equation is equivalent to 
the equation for the energy-loss time 

Ti{Ei) = r\^e{e)T{e))-^de (25) 

JO 

proposed in the papers^^i^^. 

One more function of interest is also the time of setting of the quasi- stationary electron 
distribution. Now we assume that the source of light is switched on at t=0 and slowly 
changes afterwards satisfying the condition dSo{t)/dt « So(t)B{e). The time of setting is 
determined as the time necessary for transient processes to extinct after the light is switched 
on. Employing definitions ( 1T3| |23|) and neglecting the energy dependence of P we can rewrite 
Eq. (dl]) as 

^^ = nit, e + hcoone + hu,) - n{t, e)r(e - UuJo) + ^!^!!!^^ (26) 

We solve this equation supporting on a set of excess electron energies e^ = TiUm + P^^o, P 
= -^ m, that is from hum to the maximum excess energy e^- Introducing the notation 
Tp = P{ep)N{ep) we transform Eq. fl28|) to the set of equations near the e^ levels: 

dn{t,em) I, x-p , dninst(t,em) 

-^r- = -"^'' '-^^-' + — di — 

= n(t, e^)Vm - n[t, e„_i)r^_2 H — 

at at 

On[t, €-m-p) _ fj. \-p _ fj. \-p _|_ (y'^instyt, ^m-p) 



M^ = ^(t,eOr, + ^!^i#^ (27) 

ot at 

For the highest energy level the solution is 

n{t, e„) = [ r 9n^nst{f,e^) ^^_^, ^ ^^^^ ,j^,-r^-.t (28) 



/o dt 

For the lower energy levels, except of the 0-th, the solution obtained by recursion is 

nit, e^_,) = [J^ dfi,Ut^,e^~p) ^r^.^^.t' ^ ^^q^ e^_^)]e-r--*. (29) 



Here we defined the modified source function 



+ nit, em-p+i)r.m-p+i (30) 



dt dt 

The solution for the m-th level has terms proportional to e~'"'"~i*, so omitting these terms 
decaying with characteristic time Tm = l/r~^i one obtains static solution 

n{t,em) = -Ki /^„^_l (31) 

The solution for the level em-i becomes static after the time 

Tm-l = Tm + '^/^m-2 (32) 

and it has the form 

n{t, em-i) = t; :^ [ jr. + tz: r^_i] (33) 

J- m-lJ- m-2 Ot at 

It follows from further recursion that for an arbitrary level e^ the static solution is realized 
after the time 

m—l 

rp=Y. l/r. (34) 

s=p 

and it is 

n{t, e,) = -^ E ^""-f^'-^ r^ (35) 

J- pJ- p-l s=p ot 

Taking into account the smallness of the huQ value and the definition (!23|) we can replace 
these sums with integrals 



and 



TsiE,) = (huo)-' / T{e)~'de (36) 
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Eq. (136|) is similar to Eq. (122|) . but contrary to the time of energy loss the time of EDF 
setting is determined by the time of the transitions to Ei from all the higher states. It is 
easy to show that if the source function is defined via Eq. (jl]) the Eq. (ETj) becomes, as it 
has to be, equivalent to Eq. ([7]). 

III. TECHNICAL DETAILS 

We apply the described approach to the cases of ZnO in the structure of wurtzite and 
Ti02 in the anatase structure. Pure ZnO is the semiconductor with direct band gap of 3.4 
eV, whereas anatase has the band gap 3.2 eV wide. The energy band structure of the pure 
anatase and zinc oxide has been extensively studied, mainly by the methods of the density 
functional theory. Normally such methods produce the value of band gap much less than the 
experimental data. The present study is based on the band structure calculations2Si2^i^ for 
the pure and doped anatase modified by applying single-site coulomb correlation corrections 
within the LSDA+U approach based on the LMTO band structure method. Such approach 
produces for the pure and doped anatase the values of the band gap and the energies of 
impurity states close to experimental data. However, this method of the band gap correction 
is not sufficient for zinc oxide. So in this case we apply to the conduction band states of zinc 
oxide the 'scissor operator', hence we perform a rigid shift of all the conduction band states 
to higher energy until a good value of the band gap is achieved. Such approach is justified by 
the comparison with the results of the band structure calculations of ZnO corrected by the 
application of the GW many-body theory^i. They demonstrate that the application of the 
many-body GW corrections produces almost uniform shifts of the conduction band states 
to higher energy. 

IV. RESULTS AND DISCUSSIONS 

In Figs. [1] and [2] the main data on the band states, lEDF and EDF are given for zinc 
oxide. The bands of ZnO between -6.2 and - 4 eV are those composed of the Zn 3d-states, 
the states between -4 eV and the Fermi level consist of the O 2p valence states, whereas the 
conduction band states, composed of the Zn 4s-states, are separated by the direct band gap 
3.4 eV wide. 
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FIG. 1: Left panel: the density of states for ZnO; right paneh electron dispersion curves of zinc 
oxide near the band gap. Both kinds of data are given with respect to the Fermi level. 







FIG. 2: lEDF S{e) and EDF /(e) values for zinc oxide given as functions of the excess energy of 
excited electrons e and excitation energy Eg^c ■ 

It follows from the Fig. [2]that in the case of zinc oxide the main feature of the S{u)) lEDF 
function is a gaussian peak whose maximum energy increases hnearly with the rise of the 
excitation energy. This peak corresponds to the vertical transitions from the highest valence 
band states to the lowest conduction band states for the points located at the T — K direction 
of the Brillouine zone and in the vicinity of this direction. With the increase of the excitation 
energy the wave vector of the transition shifts from the F to K point. The transition matrix 
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FIG. 3: Dependence of the relaxation rate P on the excess energy e for zinc oxide. In the inset also 
the energy dependencies of the total relaxation rate F and density of states N are given. 

500- 




FIG. 4: Dependence of the EDF setting time r^ (thick solid line) and energy loss time ti (thick 
dashed line) on the excess energy e for ZnO. Also the values of the ti (thin solid line) and Tg (thin 
dashed line) are given calculated with the constant transition probability P= 1 eV/fsec. 

elements of such excitations are much higher than the TME for excitations from all the 
lower valence band states, so the excitations from the lower states manifest themselves as 
the low satellites of the main peak. In correspondence with this lEDF function, the quasi- 
stationary function f{uj) demonstrates a Fermi-like distribution of the excited electrons in 
the conduction band, with almost constant number of electrons at the energy levels from 
the bottom of the conduction band to the energy Eexc — 3.4 eV. 
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In Figs. [3] and S] the data concerning electron dynamics in ZnO anatase are given. Fig. 
[3] also characterizes the applicability of the "effective phonon" approximation with energy- 
independent transition probability P. The change of the total probability r(e) with energy 
is almost linear beginning from the lowest excess energy e = hUm, whereas the density of 
states A^(e) also changes almost linearly, but reduces to very small values at the low limit. 
Therefore the value -P(e) = r(e)/A^(e) diverges near the bottom of the conduction band, but 
at the energy above 0.3 eV the change of -P(e) with energy is slow. 

The calculated data on the energy loss time r/ and EDF setting time r^ are shown in Fig. 
HJ They demonstrate that the relaxation of the electrons to the bottom of the conduction 
band occurs within the time less than 500 fsec. The setting of the quasi-stationary electron 
distribution in presence of the steady in time light source also occurs within the time of 
no more than 500 fsec. These data are in reasonable correspondence to the experimental 
data^~— . An analysis of the problems of comparing the experimental and theoretical data 
can be found in the Ref.— . 

In Fig. Hlalso the data on r/, r^ are given calculated with a constant value of the transition 
probability P. In order to obtain the best results the integration interval from to 0.3 eV has 
been from the calculations excluded. The calculations demonstrate that the approximation 
of the energy-independent P produces sufficiently good results only at the excess energy 
more than 0.7 eV. Variation of the P- value does not help to improve the correspondence to 
the results of the calculations with the energy-dependent probability. 

In Fig. Othe density of states and dispersion curves for anatase are given. The conduction 
band states at the energy from 3.2 to 7.8 eV above the Fermi level are composed mainly of 
3d Ti states. The corresponding density of states sharply changes with energy that evokes 
essential variation of F-probability^^. On opposite to the case of ZnO the calculations 
demonstrate that the band gap of anatase is not direct. The highest valence band state 
with almost equal energy are observed in Z- and M-point whereas the lowest conduction 
band state is in F-point. The decrease of energy of the highest valence states on direction 
from M to F points has been confirmed recently on ARPES experiments^. 

At the excitation energy from 3.5 to 3.8 eV the S{u) dependencies for anatase also are 
gaussian peaks accompanied by low-energy satellite. Similarly to the the case of zinc oxide 
they correspond to the transitions from the highest occupied valence band states to the 
lowest conduction band states. The wave vectors of these excitations are located on the 
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FIG. 5: Left panel: the density of states of anatase; right panel: energy band structure of anatase 
near the band gap. Both kinds of data are given with respect to the Fermi level. 
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FIG. 6: The dependencies of the lEDF S{e) and EDF /(e) functions on the excess energy of excited 
electrons e and excitation energy Eg^c for anatase. 

r — Z direction or in the vicinity of it. For example, the excitations at 3.5 eV occur from the 
states near the middle point of the T — Z direction, the excitations at 3.6 eV occur from the 
states near Z point, and the excitations at 3.8 eV occur from the states near the F point. 
However, contrary to the case of zinc oxide, the TME for excitations from some lower band 
states are higher than those for the excitations from the highest valence band states. Such 
are the band states that belong to the 13-th and 14,15-th degenerate dispersion curves on Fig. 
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FIG. 7: The dependence of the relaxation rate P on the excess energy e for anatase. In the inset 
also the energy dependencies of the total relaxation rate F and density of states are given. 

IHand the states of general symmetry in the vicinity of these curves. So, when the excitation 
energy increases up to the value sufficient for the excitations from these lower bands, and 
this energy is of 3.9 eV, a shift of the S{u) peaks to the lower u occurs. At energy above 
this threshold the excitations from the highest valence band states to the lowest conduction 
band states contribute to only low satellites of the the S{u) functions. In comparison with 
ZnO, this modification of S{uj) provokes essential change in the corresponding f{u) EDF 
function. At energy above 3.9 eV it is no more a Fermi-like regularity, but is composed of 
the main peak at e <0.6 eV and a low shoulder extending up to the excess energy equal to 
Eexc — 3.2 eV. So at any Eexc and e the number of excited electrons in anatase is markedly 
less than in ZnO. 

The characteristics of electron dynamics are given for anatase in Figs. [7] and [HI As in the 
case of ZnO, the transition probability P has divergence near the bottom of the conduction 
band. The interval of rapid change of P with energy extends from zero to about 0.7 eV. 
At higher energy P varies near 0.25 eV/fsec, the value much less than in the case of ZnO. 
Nevertheless, as it follows from Fig. [HI both the energy-loss time ti and EDF setting time 
Ts appear to be much less than for ZnO. This is associated with much higher values of 
density of state that provide higher rate of electron-phonon relaxation. Also the r/ and r^ 
values are shown in Fig. [8] calculated with constant value P=0.25 eV and with interval of 
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FIG. 8: The dependence of the EDF setting time Tg (thick sohd hne) and energy loss time ti (thick 
dashed hne) on the excess energy e for anatase. Also the values of the r; (thin solid line) and Ts 
(thin dashed line) are given calculated with the constant transition probability P= 0.25 eV/fsec. 

rapid variations of P excluded from integration. It is evident that these calculations only 
reproduce the trends of changes, whereas the actual deviations from the exact r/, Tg-values 
mounts to 400 %. The variations of the P- value between 0.1 and 0.4 eV/fsec and also of the 
low limit of integration do not help to improve these results. 

V. CONCLUSIONS 



We have proposed a first-principle method for evaluations of the distribution function 
of the excited electrons in the conduction band of semi-conductors. The approach takes 
into account the loss of electron energy via electron-phonon coupling, the emergence of 
new electrons excited by the external light source and redistribution of electrons between 
the energy levels. The method has been applied to the evaluations of the static electron 
distribution function in anatase and zinc oxide; the time of energy loss and the time of 
setting of the static distribution also have been calculated. 

The method helps to come also to some conclusions that may have relation to the photo- 
catalytic activity of the compounds. It is generally accepted that the photocatalytic activity 
of the oxide semiconductors is to a great extent determined by the light absorption because 
it is proportional to the number of created electron-hole pairs. Therefore the band structure 
and optical absorption in pure and doped ZnO and Ti02 was a subject of a great number 
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FIG. 9: Left panel: the energy dependencies of the /(e) and r(e) values for ZnO and anatase; 
right panel: the energy dependencies of the r(e)/(e) products. The data on ZnO are given in solid 
curves, and the data for anatase are in dash. 



of calculations, see e.g. Ref. 36|]. Also a factor important for the photocatalytic activity is 
the rate of charge transfer between the oxide and molecules absorbed on surface; this rate 
should be higher than the rate of electron-hole recombination. This factor deserves a special 
attention, but till now it was only hardly touched in the first-principle approaches. Our 
calculations permit to discuss one more characteristic of oxides that also can be important. 
Namely, one should expect that, irrespectively of the details of the processes on the surface, 
the less is the energy loss of excited electrons in bulk the more is the portion of the absorbed 
energy of light that can be spent to produce a photochemical reaction. So it is worthwhile 
to introduce the value 

ep,(e) = ^a;o/(e)r(e) (38) 

as the measure of the energy transfer to the phonon system. In Fig. [9] we compare these 
values for ZnO and anatase (the factor Tiuj^ is omitted). They are given for the excess 
energy of 0.3 eV, that is the middle energy between the band gap and the edge of the 
solar light, 4 eV. The /(w) value is in this region slightly less for anatase than for ZnO, 
but the F value for anatase is much higher. Evidently, this is associated with much higher 
density of states in the conduction band of anatase. Therefore the rate of energy loss in 
ZnO appears to be about 5 times less than in anatase. This is the factor favoring the higher 
photocatalytic activity of ZnO. Since this ratio depends on the states of only conduction 
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band and on the electron-phonon coupling of only these states, one can expect that such 
ratio is valid also for the ZnO and anatase doped with elements producing states inside the 
band gap. So we may expect that the doped ZnO can probably be more promising material 
for photocatalytic applications than the doped anatase. The photocatalytic properties of the 
doped ZnO have been studied much less than those of the doped Ti02. Nevertheless, there 
are few papers containing the comparison of activity of ZnO and anatase in the reactions of 
photo-decolarization of red dye^*^. It has been shown in^l that the ZnO catalyst composed 
of nano-particles has the activity equal to that of the analogous anatase catalyst, in spite 
of that the surface area of the anatase catalyst was more extensive. The comparison of 
the effectiveness of ZnO- and Ti02-based photo-catalysts in degradation of dyes made by 
the authors of the Ref.— demonstrated a 20 - 30 % superiority of the ZnO-based catalysts. 
There are also examples of the doped ZnO-based catalysts that have high activity in visible 
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